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ABSTRACT 

The central mass distribution of clusters of galaxies can be inferred from 
gravitationally lensed arcs with known redshifts. For the Abell cluster 2218, 
this method yields a core mass which is larger by a factor of 2.5 it 0.5 than 
the value deduced from X-ray observations, under the assumptions that the gas 
is supported by thermal pressure and that the cluster is spherical. We show 
that a non-thermal pressure support is the most plausible explanation for this 
discrepancy. Such a pressure can be naturally provided by strong turbulence 
and equipartition magnetic fields (~ 50/iG) that are tangled on small spatial 
scales 10 kpc). The turbulent and magnetic pressures do not affect the 
measured Sunyaev-Zel'dovich effect for this cluster. Intracluster magnetic fields 
with a comparable magnitude (~ 10^~^/iG) have already been detected by 
Faraday rotation in other clusters. If generic, a small-scale equipartition magnetic 
field should affect the structure of cooling fiows and must be included in X-ray 
determinations of cluster masses. 

Subject headings: gravitational lensing - galaxies: clusters of - magnetic 
fields - cosmology: dark matter 

1. Introduction 

The central mass distribution in clusters of galaxies can be inferred in an unambiguous 
way from gravitational lensed arcs with known redshifts (Grossman & Narayan 1989). 
Surprisingly, this method yields a core radius which is smaller by a factor of a few than the 
value deduced from X-ray observations (e.g., Jones & Forman 1984; Fdge & Stewart 1991) 
under the standard assumptions of thermal hydrostatic equilibrium and spherical symmetry. 
This discrepancy was highlighted recently by Miralda-Fscude & Babul (1994). Possible 
explanations to the discrepancy involve different violations of the standard assumptions, such 
as deviations from sphericity in the cluster mass distribution or non-thermal contributions to 
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the gas pressure from magnetic fields and bulk motions. In this paper we show quantitatively 
that the most plausible explanation to the above discrepancy in the rich cluster A22f 8 is the 
existence of turbulence and equipartition magnetic fields in the cluster core. 

Strong magnetic fields were detected in the past in many clusters of galaxies through 
Faraday rotation of radio sources (Kim et al. 1990, I99I; Taylor & Perley 1993). The 
detected fields have a typical magnitude of a few //G and a coherence length of I0^~^ 
kpc. Although significant, these fields are still an order of magnitude smaller than the 
equipartition values needed to make them dynamically important for the support of the 
cluster gas. However, the Faraday rotation method is limited from probing fields that 
are tangled on small spatial scales, because field reversals along the line of sight cancel 
out in the observed rotation measure. The possibility therefore remains that somewhat 
stronger magnetic fields, tangled on small scales 10 kpc), make a significant contribution 
to the pressure support of the cluster gas. Indeed, evidence for tangled magnetic fields with 
magnitudes as high as ~ 10^~^/iG was found in several clusters (Dreher et al. 1987; Perley 
& Taylor 1991; Taylor & Perley 1993; Ge & Owen 1993). Such fields should result in a 
discrepancy between the thermal pressure of the gas and the depth of the cluster potential 
well. In this work we argue that this discrepancy was in fact detected through gravitational 
lensing in the Abell cluster 2218. 

The outline of this paper is as follows. In §2 we formulate the hydrostatic equilibrium 
equations for a hot gas with an isotropic (small-scale) magnetic stress, embedded in a cluster 
with a general ellipsoidal geometry. We then infer empirically the density distribution of the 
gas according to the observed profile of its bremsstrahlung surface brightness. At hydrostatic 
equilibrium, the fitted gas distribution can then be used to find the net gravitational potential 
of the cluster. Gravitationally lensed arcs probe the total cluster mass enclosed within 
a cylinder at some projected radius. This mass can be related to the observed surface 
brightness of the X-ray gas using Poisson's equation. This empirical approach (cf. Fabricant, 
Rybicki & Gorenstein 1984) has the attractive feature that it does not enforce any particular 
form of the total mass distribution; instead, the lensing mass is obtained directly from 
the observed X-ray surface brightness profile. The method is applied in §3 to the Abell 
cluster 2218 for which there are a number of additional observational data, including the 
velocity dispersion of the member galaxies (Le Borgne et al. 1992) and measurements of the 
Sunyaev-Zel'dovich (SZ) effect (Birkinshaw & Hughes 1994). We show that the mass of the 
projected cluster core has to be larger by a factor of 2.5 it 0.5 than the value inferred from the 
standard X-ray analysis. We then demonstrate that this discrepancy cannot be explained 
by an elongation of the cluster along the line of sight. On the other hand, the discrepancy 
can naturally result from strong turbulence and equipartition magnetic fields in the cluster 
core. In §4 we summarize these conclusions and discuss ways to test them empirically. 
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2. Empirical Relation Between The Lensing Mass and X-Ray Data 



The standard description of the hot gas in clusters of galaxies assumes hydrostatic 
support by thermal pressure and spherical symmetry. To consider deviations from these 
standard assumptions we express the total gas pressure as the sum of thermal and non- 
thermal components p = pt + Pnt^ and assume a general ellipsoidal geometry for the cluster 
gas. The thermal pressure of the gas at a temperature T is given by pt = PgkT / finip^ where 
Pg is its mass density, k is the Boltzmann constant, nip is the proton mass, and p 0.6 
is its mean atomic weight. The non-thermal pressure may contain three separate parts 
(which are comparable in magnitude for the interstellar medium of the Galaxy), namely: 
magnetic pressure p^, turbulent pressure pturb^ and cosmic-ray pressure Pray A magnetic 
field that is tangled on scales small compared to the cluster core radius (~ 100 kpc) yields 
a diagonal magnetic stress tensor with a pressure^ ps = {\B\'^) /Sir. Turbulent motions and 
bulk velocities add pturb = \{Pg'v'^)^ while cosmic-rays confined by the magnetic field provide 
their own kinetic pressure p^ay ^ Pb- 

The hydrostatic equilibrium equation is 

^Vp=-V$, (1) 

pg 

where $ is the gravitational potential of the cluster. Since gravitational lensing probes 
mainly the cluster core we assume for simplicity that 

Pnt = OL Pt, a = const; (2) 

and that the gas is isothermal. Equations (1) and (2) yield. 

The total mass density of the cluster can then be inferred from Poisson's equation, 

1 „o (1 -|- a)A;r„„ , -I , , 

Pioi = = ^ ^ ' V In p;' , 4 

47rG AirGpnip ^ 

if the gas density distribution is known from X-ray observations. 



"'^The magnetic pressure contains two components, arising from the magnetic energy density Mmag = 
(B^/Stt) and the kinetic energy density t/kin = ^mag of Alfven waves (Dewar 1970). Other modes 
are Landau damped (Kulsrud 1994). The magnetic virial theorem (Shu 1992) yields the total pressure 

Pb = fukm + ^Mmag = / StT . 
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The X-ray surface brightness of the cluster gas is obtained by a hne integral of its 
bremsstrahlung emissivity, 

I oc j pI dz, (5) 

where we use cylindrical (r, z) coordinates around the line of sight. Typically, the observed X- 
ray brightness contours are close to circular. We therefore assume axial symmetry around the 
line of sight, and adopt the simplest model to describe deviations of clusters from sphericity. 
In particular, the gas density is assumed to be a function of the ellipsoidal coordinate, 

„2 ^2 



,2 



r z 



m 



(6) 



r1 z'^ 



where (tc, Zc) is the ellipsoidal core size. For the simple case of a spherical gas distribution, 
m = R/Rc is the radial coordinate and Rc = Vc = z^ is the core radius. With pg = Pg{m) we 
can use Abell's equation (e.g., Binney & Tremaine 1987) to invert equation (5), 



2 ( dl\ dr 



Typically, the X-ray surface brightness profile /(r) is well fitted by a functional form (e.g. 
Birkinshaw, & Hughes 1994; Bahcall & Lubin 1994, and references therein) 

^Mo.(l + ;j) - (8) 

with /3 = 0.5 — 0.9. Equation (7) then yields a gas distribution of the form, 

/>,(m) oc (l + m2)-3'^/2. (9) 



Gravitational lensing probes the total mass enclosed within a cylinder of a particular 
radius r = b along the line of sight. To find this mass we calculate the volume integral, 

/ , (l + a)kT / - , (l + a)kTb , dlnp„ , , 

where we have substituted equation (4), and used Gauss's theorem to convert the volume 
integral into a surface integral on an infinite cylinder. Note that the lensing mass 
Af(6) depends only on the shape of the gas density distribution and not on its absolute 
normalization. Using equation (9) we finally obtain 

'nhT 38 P , , Zr , , 

M(b) = - , (1 + a) —. 11 
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This is our basic result. Observational measurements of T, /3, rc, b and M[b) can be used to 
determine the quantity (1 + a)zc/rc = rj. Values of rj ^ 1 quantify the degree to which the 
standard assumptions concerning the cluster g violated. 

Finally, we consider another probe of the thermal state of the cluster gas, namely 
the Sunyaev-Zel'dovich (SZ) effect. This effect results from Thomson scattering of the 
cosmic background photons by hot electrons in the optically-thin cluster. If the cluster 
has no peculiar motion, the single-scattering Doppler shift averages to zero to linear order 
in the electron velocity, since the distribution of electron velocities is isotropic at thermal 
equilibrium. However, the second-order Doppler effect does not vanish and is proportional 
to the square of the thermal velocity of the electrons, i.e. to kT. The scattering probability 
is proportional to the integral of the electron density along the line of sight, and thus the net 
SZ effect is linear in the electron pressure. The scattering of microwave photons to higher 
energies leads to a decrease in the Rayleigh- Jeans temperature of the cosmic background 
radiation Trj, ATrj/Trj oc — / pfdz. Subsonic turbulent velocities provide vf^^^^ ^ (kT/nip) 
and result in an effect that is smaller than the thermal effect of the electrons by the electron 
to proton mass ratio. The existence of magnetic fields and cosmic-rays can also be ignored 
in calculating the SZ distortion. Thus, for the density profile of equation (9) the distortion 
is given by, 

/ 2 \ 1/2-3/3/2 

ATnj{r) « (^1 + ^ ^, (12) 

with no reference to the non-thermal pressure of the gas. 

Note that according to equations (8) and (12), the elongation of clusters does not change 
the radial shapes of either the X-ray surface brightness or the SZ distortion. The elongation 
only enhances both effects by the axis ratio factor Zc/vc- 

3. Application to A2218 

Next we apply the results from §2 to study the conditions of the gas in the rich cluster 
A22I8, for which there is a wealth of data, including the X-ray surface brightness profile of 
the cluster, the positions and redshifts of gravitationally lensed arcs, the Sunyaev-Zel'dovich 
effect, and the velocity dispersion of the member galaxies. We denote the Hubble constant 
by Ho = 100 h km s'^ Mpc-\ 

Abell 2218 is a cluster at a redshift z = 0.175 (Le Borgne et al. 1992) with a richness 
class 4, an X-ray luminosity in the 2-IOkeV band of 10"*^ erg s~^ (David et al. 1993), a 
central proton density of (5.4 it 0.5) X 10~^[h/0.5y^^[zc/rc)~^^^cni~^ (Birkinshaw & Hughes 
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1994), and an X-ray temperature of kT = 6.7to!4 keV (McHardy et al. 1990). At the 
redshift of the cluster, 1 = 1.9 h kpc. The optical image of this cluster reveals several 
gravitationally lensed arcs. One of the arcs is at a redshift of Zs = 0.702 with a critical 
angular radius of 9crit = 20.8 centered at the X-ray peak of the cluster (Pello et al. 1992; 
see also Miralda-Edcude & Babul 1994). The SZ effect for A2218 was most recently analyzed 
by Birkinshaw & Hughes (1994), who improved previous work by McHardy et al. (1990) 
and derived a Hubble constant of h = (0.65 it 0.25) X [rd Zc). Birkinshaw & Hughes (1994) 
fitted the X-ray surface brightness profile using equation (8) and obtained the best fit with 
/3 = 0.65 lb 0.04, and an angular core radius of = 60 it 10 . With these values, equation 
(8) also provides an excellent fit to the surface brightness profile for this cluster as recently 
measured by ROSAT (Stewart et al. 1993). We therefore adopt these values in the analysis 
that follows. 

The observed arcs primarily constrain the total mass inside the cylinder of the critical 
radius h. The constraint is relatively insensitive to the mass distribution profile (e.g., 
Kochanek 1991), and can be expressed as 

where D^^Ds, and D^s are the angular diameter distances to the lens, to the source, and 
from the lens to the source, respectively (see Kochanek 1992 for the redshift dependence of 
D on different cosmologies). Note that h = OcritDL and Vc = OcDl- For A2218, equation (13) 
gives M[b) = 0.32h~^ x lO^'^Af© assuming a cosmological density parameter = 1. This 
result has a negligible dependence on the underlying cosmology. In fact, for any cosmology 
with < < 1 and no cosmological constant (A = 0), or for any fiat cosmology with A 7^ 0, 
this mass constraint varies by less than 3%. Equation (13) assumes an axially-symmetric 
lens geometry. In reality, the contribution from the presence of a central cD galaxy and 
substructure in A2218, and the ellipticity needed to break the cylindrical symmetry and 
form arcs instead of rings, lower this mass estimate by 15% (Miralda-Escude 1994). 

However, from equation (11) we obtain M(6) = (0.11±0.02)(A;r/6.7keV)?7 h'^xlO^^MQ. 
In order for these two mass estimates to be consistent, we find 

?7 = (1 + a)— = 2.5±0.5, (14) 
r-c 

where the error bars are dominated by the uncertainties in the X-ray observations and can 
be reduced considerably in the future. 



The first possible origin for this value of rj is that the cluster is prolate. If the contribution 
of non-thermal pressure is negligible, then the axis ratio for the gas distribution must be 
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Zc/vc ~ 2.5. We then find, using equations (4) and (9), tliat tlie total mass distribution has 
an axis ratio of roughly 6:1 in the cluster core. The axis ratio for the total mass is larger than 
Zc/vc because the gas follows the potential which tends to be smoother than the underlying 
mass distribution (Binney & Tremaine 1987). Furthermore, the total mass distribution has 
an unphysical dumbbell shape at large radii for any value of [zc/rc] > y^3/2, as long as 
the gas distribution is described by equation (9) (cf. Kassiola & Kovner 1993). The above 
properties of the total mass distribution do not seem to be plausible. With an axis ratio 
6:1, the probability for a perfect alignment of the cluster along the line of sight (as required 
by the circular X-ray brightness contours) is only a few percent times the relatively small 
fraction of prolate clusters seen from the side (Jones & Forman 1992). In addition, such a 
configuration is unstable against bending modes (Merritt & Hernquist 1991). Therefore, a 
full account for the discrepancy by the prolateness of the cluster is unlikely. 

The alternative explanation involves a non-thermal pressure support of the gas. This 
explanation is consistent with other observational data. The observed velocity dispersion of 
galaxies (Le Borgne et al. 1992) 1370ti2okms-i is larger than {kT/fimpY/^ = 1050t3okms-i. 
In addition, the most recent determination of the Hubble constant from the SZ effect in 
A2118 (Birkinshaw & Hughes 1994) did not yield a value smaller by a factor of 2.5 from the 
range 0.5 ^ h ^ 1^ as would be expected according to equation (12) if the gas distribution 
was elongated along the line of sight rather than being supported by a non-thermal pressure. 

In principle, all the potential sources of non-thermal pressure are equally viable in 
accounting for the lensing mass discrepancy. However, large bulk velocities can be excluded 
by other considerations. Supersonic rotational velocities would have fiattened the X-ray 
brightness contours, while large radial velocities would have relaxed to equilibrium within a 
dynamical time (^ 10^°yr). The existence of a strong turbulent pressure pturb ~ Pt is likely 
to be accompanied by equipartition magnetic fields, in analogy with the conditions in the 
interstellar medium of the Galaxy. We therefore conclude that equation (14) implies the 
existence of dynamically important turbulence and magnetic fields in A2218. 

For example, a 2 may correspond to ps ~ Pturb ~ Pt- This amounts to an 
equipartition magnetic field strength 

« = (sTTFwi) (Tk^j • 

where Ue is the electron density. If this field is tangled on relatively small scales (^ lOkpc), it 
would partially avoid detection by Faraday rotation measurements. Previous observations of 
other clusters (Dreher et al. 1987; Kim et al. 1990, 1991; Perley & Taylor 1991; Taylor 
& Perley 1993; Ge & Owen 1993) found evidence for fields that are somewhat weaker 
or comparable to the equipartition value (10°~^/iG) on large scales (10°'^~^ kpc). Fven 
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relatively weak fields should eventually approach equipartition values at some radius as they 
are dragged inwards by the cooling flow of the cluster (Soker & Sarazin f990). 

The existence of strong fields results in synchrotron cooling of the cluster electrons. 
However, since the plasma frequency of the cluster, fp = 0.64 kHz (ne/5 X f 0~"^cm~"^)^/^, 
is larger than the cyclotron frequency, fc = 0.f4 kHz (i?/50/iG), the synchrotron emission 
is suppressed. The plasma cutoff is effective across the entire region where the field is in 
equipartition with the thermal pressure of the gas, because the ratio between the above 
frequencies scales as p^^^/B oc {pti PbY^^ ■ The synchrotron cooling time is therefore 
much longer than the bremsstrahulg cooling time in the cluster core, and any residual 
emission would occur at the unobservable frequency regime ~ kHz. On the other hand, 
the synchrotron emission by cosmic rays is observable at microwave frequencies. From the 
fact that the measurement of the SZ effect (Birkinshaw & Hughes f 994) was not dominated 
by a synchrotron signal we conclude that p^ay ^ Pt- This leaves only magnetic fields and 
turbulence as viable sources of non-thermal pressure in A22f 8. 

4. Conclusions 

The core mass of the cluster of galaxies A22f8 can be probed by two independent 
methods: gravitational lensing and X-ray observations. Assuming that the cluster is spherical 
and that the hot gas is supported by a thermal pressure, these two methods provide values 
that are different by a factor of 2.5 it 0.5. This conclusion is in agreement with a recent study 
by Miralda-Escude & Babul (f993). The result is not sensitive to the underlying cosmology. 
The main uncertainty is associated with the temperature at the cluster core, which was 
not resolved by early X-ray spectroscopy observations (McHardy et al. f990). However, a 
recent observation by ROSAT shows a fiat temperature profile, consistent with the modeling 
assumption of an isothermal cluster (Stewart & Edge f993). The above mass discrepancy 
can be resolved by relaxing either of the two underlying assumptions about A22f 8. In §3 we 
showed that in order to explain the discrepancy by a deviation of the cluster from sphericity, 
its mass distribution must be elongated along the line of sight with an axis ratio of 6:1. Such 
a distribution is highly improbable, and also gravitationally unstable (Merritt & Hernquist 
I99I). 

The alternative explanation involves a non-thermal pressure support of the gas. Large 
rotational velocities are excluded based on the circular X-ray brightness contours. Supersonic 
radial velocities would have been damped out within a dynamical time ~ 10®~^yr. Indeed, 
the X-ray map of the cluster core does not show large scale inhomogeneities of the gas. 
In addition, cooling fiow velocities are expected to be highly subsonic at the core radius. 
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However, there are other more plausible sources of non-thermal pressure. Turbulence and 
magnetic fields in equipartition with the thermal energy of the gas can fully account for 
the above discrepancy in a spherical cluster. This explanation is consistent with other 
observational data on A2218. In difference from a cluster elongation, the turbulent and 
magnetic pressures do not affect the Sunyaev-Zel'dovich effect, in line with the observations 
(Birkinshaw & Hughes 1994). The square of the observed velocity dispersion of galaxies (Le 
Borgne et al. 1992) is larger than {kT/finip) by a factor of I.Ttois! a value that may also 
refiect the different bias and virialization history of the galaxies relative to the gas in the 
cluster. 

Magnetic fields with an amplitude ^ fiG were detected by Faraday rotation in many 
other clusters (Dreher et al. 1987; Kim et al. 1990, I99I; Owen et al. 1990; Perley & Taylor 
I99I; Taylor & Perley 1993; Ge & Owen 1993). In the Hydra A cluster, Taylor & Perley 
(1993) found a magnetic field ~ 6/iG on a scale ~ 100 kpc, and a tangled field as high 
as ~ 30/iG on smaller scales. In the cluster AI795, Ge & Owen (1993) observed rotation 
measures exceeding 3000 rad m~^, which translate to a field amplitude > 20/iG. Similar 
results were obtained by Dreher et al. (1987) for the Cygnus A cluster, and by Perley & 
Taylor (I99I) for the 3C295 cluster. The tangled magnetic field we predict for A22I8 is 
comparable in amplitude to the inferred fields in these cases. Although the average rotation 
measure (RM) of such a tangled field is close to zero, the dispersion in the rotation measure 
should be large, (RM^)!/^ = 2I40(ne/5 x I0-3)(5/50/iG)(£/I0kpc)(iV/I0)^/2 ^^d m'^, where 
N is the number of cells of size £ and field strength B along the line of sight. 

Aside from the existence of strong magnetic fields, our model predicts a considerable 
discrepancy" (cf. Bahcall & Lubin 1994, and references therein) for A22I8. The radial 

distributions of the gas and the galaxies should not be consistent with their relative velocity 

dispersions due to the non-thermal support of the gas. 

Although galaxy-driven turbulence cannot readily produce the observed fiG magnetic 
field (De Young 1992), there are a variety of field generating mechanisms operating in the 
cores of clusters (Taylor & Perley 1993). A natural source for the needed fields is dynamo 
amplification of seed galactic fields (Ruzmaikin et al. 1989; Kulsrud & Anderson 1992) due 
to the strong turbulence and shearing present during the collapse and virialization stages 
of the cluster. It is interesting to note that in an = I universe, the accretion of clumps 
of matter by the cluster continues at all times (Richstone, Loeb & Turner 1992) and can 
persistently excite turbulence in it. In analogy with the interstellar gas, the amplification 
of the magnetic fields would saturate at equipartition with the turbulent pressure. The 
resulting fields would then affect the structure of the cooling fiow in the cluster (Soker & 
Sarazin 1990). 
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If future observations show that turbulence and magnetic fields are dynamically 
important not only in the inner parts but also in the outer parts of rich clusters, this would 
have a variety of interesting implications. First, it would increase the cluster mass to light 
ratio by a factor of a few and shift the resulting estimate of from the range 0.2-0.3 (Binney 
& Tremaine 1987) to about unity. The baryonic mass fraction of clusters would consistently 
decrease to values ^0.1. Both of these changes would tend to make the cluster data agree 
with the predictions of infiation and standard nucleosynthesis. Second, an increase in the 
high mass tail of the cluster mass function would affect the comparison between observational 
data and popular cosmological models (Bahcall & Cen 1993). Finally, this would raise the 
possibility that magnetic fields have a non-negligible infiuence on structure formation in the 
universe. 

We thank Daniel Fisenstein, Alyssa Goodman, Russell Kulsrud, Jordi Miralda-Fscude, 
Richard Mushotzky, and Ramesh Narayan for useful discussions. 
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